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Abstract 

The pure leptonic decay fi^ ji~ is strongly suppressed in the Stan- 

dard Model (SM), but can have large enhancements in Supersymmetry, 
especially at large values o/tan/3. New limits on this decay channel from 
recent LHC data have been used to claim that these limits restrict the 
SUSY parameter space even more than the direct searches. However, 
direct searches are hardly dependent on tan/3^ while Br{B^ l^^ 
is proportional to tan^/3. The relic density constraint requires large 
tan P in a large region of the parameter space, which can lead to large 
values of Bg ji^ ii~ . Nevertheless, the experimental upper limit on 
Br(Bg is not constraining the parameter space of the CMSSM 

more than the direct searches and the present Higgs limits, if combined 
with the relic density. We also observe SUSY parameter regions with 
negative interferences, where the B^ jji^ iJi~ value is up to a factor 
three below the SM expectation, even at large values o/tan/3. 

1 Introduction 

Flavour Changing Neutral Currents (FCNC), like the leptonic decays of neutral B-mesons, are strongly 
suppressed in the Standard Model (SM), since they can only occur via loops involving the weak bosons. 
These decays are helicity suppressed, so the amplitudes are proportional to the mass of final state 
particles and the highest rates will be into tau leptons. The experimental signature for leptonic decays 
is clear: search for an invariant mass in the mass window of the B-meson. This is easier for muonic 
decays. Hence, muonic B-decays have been investigated in much more detail at hadron colliders, 
especially since these decays can be strongly enhanced by loop corrections involving particles beyond 
the SM, like Supersymmetry [LiSj. The B^ 11^ ii~ decay mode has received significant attention 
[ZHS] after the CDF collaboration announced a measurement a factor five to six above the expected 
SM value [lOj. However, the excess was not confirmed by subsequent LHC measurements [Uj, but 
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Figure 1: Annihilation diagrams for the hghtest neutrahno, which is a hnear combination of the 
gaugino and Higgsino states: \xo >= A^il^o > ^N2\W^ > ^NslHi > -\-N4\H2 >. The 
dependence of the amphtudes on masses and neutrahno mixing parameter Ni has been 
indicated. 



nevertheless the LHC upper limit can give significant constraints on the SUSY parameter space, see 
e.g. jl2| . where in some scenarios better limits than those obtained from direct searches have been 
claimed. However, the excluded parameter space depends strongly on the choice of tan /3 since the 
rate varies as tan^/3. The relic density constraint correlates tan^S with the SUSY mass 
parameters [13], so if one combines the cosmological constraint with the accelerator constraints there 
is no arbitrary choice for tan /3 anymore. Although the relic density requires a large value of tan /3 in 
a large region of parameter space, we show that the excluded SUSY mass ranges are well below the 
LHC constraints from direct searches [TH - FTT] and the limits on the pseudo-scalar Higgs [HI [19]. In 
principle, other constraints, like g-2, h ^ and B ^ rv could also be considered. However taking 
these into account requires a careful treatment of the non-gaussian systematic errors, which is beyond 
the scope of the present letter. Numerous studies combining these variables with the recent LHC data 
have appeared [20r.22j . 

2 Relic Density 

The relic density and annihilation cross section a are related through: 

nu^ = (1) 

<(JV> ^ ^ 

where the annihilation cross section a averaged over the relative velocities of the neutralinos is given in 
pb [SSlEl] and h ^ 0.71 is the Hubble constant in units of 100 {km/ s) /Mpc. The best value for the relic 
density is ^h? — 0.1131±0.0034 [25j. For a given relic density the annihilation cross section is known 
independent of a specific model, since it only depends on the observed Hubble constant and the ob- 
served relic density. Its value is furthermore largely independent of the neutrahno mass (except for 
logarithmic corrections) [231 [21] . The DM constraint should exist for any model, but to be specific the 
Constrained Minimal Supersymmetric Standard Model (CMSSM) with supergravity inspired breaking 
terms, will be considered [26H28]. It is characterized by 5 parameters: mo, ^1/2? tan/3, sign(/x), Aq. 
Here mo and mi/2 are the common masses for the gauginos and scalars at the GUT scale, which 
is determined by the unification of the gauge couplings at this scale. Gauge unification is perfectly 
possible with the latest measured couplings at LEP [29j. Electroweak symmetry breaking (EWSB) 
fixes the scale of ji [30], so only its sign is a free parameter. The positive sign is taken, as suggested 
by the small deviation of the SM prediction from the muon anomalous moment. 

The relic density can be calculated from the diagrams in Fig. [T] For its calculation we used the 
public code micrOMEGAs 2.4 [31] [32] combined with Suspect 2.41 as mass spectrum calculator 
The optimal parameters were found by minimizing the function using the Minuit program 
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For heavy SUSY masses the sfermion exchange diagram is suppressed, the W- and Z-final states 
from t-channel chargino and neutrahno exchange have a smah cross section, the couphng of the LSP to 
the Z-boson is only via the Higgs component of the LSP, which is typicahy smah, so in most regions of 
parameter space the pseudo-scalar Higgs exchange is dominant, except for the co-annihilation regions. 
These are the regions, where the Next-to-Lightest Supersymmetric Particle (NLSP) and LSP are 
nearly mass-degenerate. In this case they co-exist in the early universe until the common freeze-out 
temperature and can co- annihilate. This happens if the stau and neutralino are degenerate and co- 
annihilate into a tau ^35j. For this to happen the values of mo and mi/2 have to be fine-tuned to a high 
degree, so it happens only in a thin stripe in the mg — ^1/2 plane, as will be shown below. Another co- 
annihilation region happens at the border of parameter space, where electroweak symmetry breaking 
does not occur anymore, since here the Higgs mixing parameter becomes negative. In the transition 
region becomes small and the lightest chargino and lightest neutralino become nearly degenerate 
Higgsinos, as is obvious from the mass matrices, which have as lowest eigenvalues either a gaugino 
mass term or a Higgsino mass term, if the mixing is neglected. In this case gauginos can co-annihilate 
into a W-boson [36j . 

Outside the bulk region with low SUSY masses and the co-annihilation regions the dominant 
contribution comes from A-boson exchange: x + X ^ ^ ^ which is proportional to 

m^ml tan^ /3 ^Nsi sin /3 - 7V41 cos {N21 cos 9w - Nu sin 9w)^ 
< av >^ 1 K ^ . (2) 

The elements of the mixing matrix in the neutralino sector define the content of the lightest neutralino: 

IX?) = Nii\Bo) + N2i\Wi) + N:,i\Hi) + Nii\H2). 

The sum of the diagrams should yield < crv >— 2 • 10~^^ crv? I s to get the correct relic density, which 
implies that the annihilation cross section a is of the order of a few pb. Such a high cross section 
can be obtained only close to the resonance, i.e. vtla ^ 2m^. Actually on the resonance the cross 
section is too high, so one needs to be in the tail of the resonance, i.e. vtla ^ 2.2m^ or rriA ^ 1.8m^. 
So one expects rriA oc mi/2 from the relic density constraint. This is shown in the left panel of Fig. 
|2j Here we optimized simply tan/3 for each pair of mg — mi/2 values, as was done in Ref. [13j. The 
corresponding values of tan /3 needed are shown in the right panel of Fig. [2j The production cross 
section of the pseudo-scalar Higgs at the LHC is proportional to tan^ /3, so the present limits from 
LHC are proportional to tan/3^ as well. At tan/3=50 the present limit on rriA is about 450 GeV, so 
the limit on mi/2 is close to it. The exclusion on ruA as function of tan /3 from the CMS collaboration 
[19] is indicated in Fig. [2] as well. Similar limits have been obtained by the ATLAS collaboration [18]. 
Also the excluded region from the direct searches has been indicated using the CMS data [16]. Similar 
results were obtained by other searches ^UlinilT^. The relic density constraint can be fulfilled with 
the parameters of Fig. [2j as demonstrated in the left panel of Fig. [3| The top left is excluded, since 
here the LSP is not a neutral particle, but the stau is the LSP. In the co-annihilation regions the 
annihilation via the pseudo-scalar Higgs exchange has to be suppressed, thus requiring a larger value 
of rriA and a corresponding lower value of tan/3, as demonstrated in the right panel of Fig. [3j It is 
just a more detailed plot of the right hand panel of Fig. [2] for two values of mi/2. 

In summary, if one allows tan/3 to vary in the mo — Tni/2 plane, one obtains the observed relic 
density for any combination of mo and mi/2, i-e. the relic density allows all masses for the SUSY 
sparticles. However, the fi^ ii~ constraint has to be investigated for the large values of tan /3 

required by the relic density. 
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Figure 2: The value of the pseudo-scalar Higgs mass (left) and the value of tan/3 required for a 
correct relic density in the mo — mi 12 plane. The excluded region by the limit on the 
pseudoscalar Higgs mass from Ref. [19j is indicated by the white solid line in the left 
panel, while the dashed line indicates the limit from the direct searches from Ref. |16| . 
The tan /3 value is around 50 in the central region, as indicated by the colour coding in the 
right panel and decreases towards the edges, where co-annihilation starts to be important. 
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Figure 3: The relic density in the mo — mi/2 plane (left) and tan /3 as function mo for different values 
of mi/2 (right) The colour coding in the left panel shows that the relic density constraint 
can be fulfilled for all SUSY masses. 
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3 ^ fi^fi decay rate 



The branching ratio for ii^ ii is taken from Ref. [6], which we write in the form 
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where fs^ is the Bg decay constant, is the B meson mass, tb is the mean hfe and mi is the mass 
of lepton. Ca, are largely determined by the SM diagrams, while Cs^ C^, Cp, Cp include the 
SUSY loop contributions due to diagrams involving particles such as stop, chargino, sneutrino, Higgs 
etc.. For large tan/3 values the dominant contribution to Cs can be written as: 



Cs 



VtbVt 



ts 



tan^/3 
4 sin^ 0w 



/ 1711)171 ij^mtfi\ sin 20^ 
M2 Ml ) 



171- 



\ 



I log 


r 2 1 
m~ 


m~ loe: 




m2 _ rjjj2 2 _ ^2 

ti 12 



(4) 



where m^^ ^ are the two stop masses, and % is the rotation angle to diagonalize the stop mass matrix. 



We need to multiply the above expression by 1/(1 + 65)^ to include the SUSY QCD corrections, where 
65 is proportional to /xtan^S [37j. We have Cp = —Cs^C'^ = {ms/mi,)Cs and Cp = —{ms/mi,)Cp. 
One observes from Eq. |4]the tan^/3 dependence, but one also observes the strong suppression in the 
last term if the stop masses become equal. In the MSSM the stop mass splitting is given by (see e.g. 
reviews 
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where the left- and right-handed quark masses are defined by: 
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rn'q + mi + -{4M^ - M^) cos 2/3, 
ml + mi - ^{M^ - Ml) cos 2/3. 



For large SUSY scales the mass terms for the right-handed singlet mu and left-handed doublet mg 
become large and m^L and mtR become of the same order of magnitude. Then the stop splitting 
is determined by the term At — ji/ tan/3, so for large tan/3 the second term is small and the stop 
mixing can be made small by increasing the trilinear coupling at the GUT scale. One indeed 
can eliminate the tension between the large value of tan/3 required by Q.h^ and the B^ 
rate, as demonstrated in Fig. [5} in the left (right) panel the dependence of Br{Bg 
Q.h? are shown as function of tan/3 for A^ — [A^ > 0). The left and right vertical scales are for 
Br{B^ IJi^lJi~) and Q.h?^ respectively and the scales have been adjusted so, that the horizontal 
line indicates the upper limit for Br(Bg ^he observed value for Qh^. One observes 

from the left panel that for the correct value of tan/3=50 for Qh'^ the value of Br{B^ 1^'^ is 
far above the experimental upper limit, but if one adjusts A^ both can be brought into agreement 
(right panel). Here we fitted simply and tan/3 for each value of mg and mi/2 i^i the mg — ^1/2 
plane with B^ ji^ ^h? as constraint. The fitted values of A^ reduce the stop mass to 
low enough values to force agreement. The required values of and the corresponding stop mass 
differences are shown in Fig. [6) The excluded regions in the combined fit of the relic density and 
Br{Bg 11^ are shown in Fig. [t] for the present limit (left panel) and for a hypothetical limit 
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Figure 5: The tan^ dependence of Br{B^ and the rehc density for Aq^O (left) and Aq > 

(right). The left and right vertical scales are for Br{Bg M^M~) and fi/i^, respectively 
and the scales have been adjusted so, that the horizontal line indicates the upper limit for 
Br{B^ fJi^ lJi~) and the observed value for Q.h?. 
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Figure 6: The Br{Bg M^M~) constraint can lead to tension in combination with the relic den- 
sity constraint, since the latter requires large tan/3, which leads to a large stop splitting. 
However, this can be compensated with a large value of Aq (left panel), which reduces 
the difference between the stop masses At (right panel) in the region where otherwise the 
constraint Br{B^ < 4.7 • 10"^ could not be fulfilled. 
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Figure 7: Excluded region from a combined fit of the relic density and the upper limit on Br{B^ 

li-^lJi-) < 1.1 • 10"^ (left) and a hypothetical Br{B^^ l^'^ < 0.66 • 10"^ (right). 
The colour code indicates the value. x^=5.99 indicates the 95% C.L. contour, which 
essentially corresponds to the red region. 
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Figure 8: Br{B^ l^'^ and the pseudoscalar Higgs mass ttia as function of and tan/3 in 
the left and right panel, respectively. The figures corresponds to mo = 1000 GeV and 
mi/2 = 250 GeV, which is inside the excluded region on the right hand side of Fig. [7j 
Note that the green region in the left panel corresponds to the SM value, while the blue 
(red) region corresponds to values below (above) the SM value. 
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of twice the SM value (right panel). One observes that the limit from Br{B^ 1^'^ is well below 
the limits from the direct and Higgs searches shown in Fig. [2j The reason for the two-lobed excluded 
regions is the following: at small values of mo the trilinear coupling cannot be made large enough to 
suppress Br{B^ l^^ enough, because the staus become tachyonic. At intermediate values of mo 
the trilinear couplings can be made large enough, but at larger values of mo the pseudoscalar Higgs 
boson mass ttia becomes too large for large values (see Fig. [s] right) and the relic density becomes 
too large as well. For values of mo well above 1 TeV the loop contributions are suppressed enough 
to fulfill the Br{Bg M^M~) constraint. These results are demonstrated in Fig. [sj which displays 
the values of Br{B^ l^^ and the pseudoscalar Higgs mass in the Ao-tan/3 plane for mo = 1000 
and mi/2 = 250 GeV, i.e. in the excluded lobe on the right hand side in Fig. [71 The green region 
in the left pannel of Fig. [s] corresponds to values close to the SM value for Br\Bg but 
at large positive values of and large values of tan/3 the Br{B^ /^^M~) value drops below the 
SM value (blue upper right region), while at lower values of Aq one observes the famous large tan/3 
enhancement (red bottom right region). In the right top corner the staus become tachyonic, so this 
theoretically disfavored region is left white. Surprisingly, values of Br{B^ l^^ can fall up to a 
factor three helow the SM value, which can be explained as follows. In Eq. |4]sin(20^) can change sign, 
depending on the value of the off-diagonal element in the stop mixing matrix At — \ij tan/3. Hence, 
Cp can change sign as well and the term \{Cp — Cp)/{mi) + m^) + 2mi^/m\^{CA — in Eq. 3 

can become small, if Cp and Ca have opposite sign. We have checked that this change in sign is 
indeed the origin of the negative interference between the SM value and the SUSY values, both in the 
micrOMEGAs code, which we used, and in the Superlso V3.1 code [40j, which gives almost identical 
results. 

4 Summary 

We have calculated the excluded regions in the CMSSM from the recent upper limits on the B^ 
ji^ ji~ decays in combination with the relic density constraint. The latter requires large tan/3 values 
in the regions outside the co-annihilation regions and since Br{B^ l^^ is proportional to tan^ /3 
one could expect strong constraints from the recent upper limits. However, the Br{Bg M^M~) 
approaches zero in case the splitting between the stopl and stop2 masses approaches zero. This 
splitting is determined by the off-diagonal element At — fi/ tan/3 of the stop mixing matrix, which 
can be made small for large tan /3 and a positive value of the trilinear coupling Ao at the GUT scale. 
From a simultaneous fit of Ao and tan/3 to the combined data of Br{B^ l^^ ^ind relic density 
we find the excluded regions from these constraints to be well below the constraints from the Higgs 
searches and direct searches at the LHC. This holds even in the case that a hypothetical limit on 
Br{Bg /i^M~) of two times the SM value would be obtained. 

It is also shown that at large values of both, tan /3 and the trilinear coupling, negative interferences 
can lead to Br{Bg M^M~) values a factor three below the SM value, so even if values below the SM 
are found experimentally, this does not exclude supersymmetry, but constrains the parameter space. 
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